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Abstract Pscudorapidity gap distributions in proton- 
proton collisions at y/s = 7 TeV are studied using a 
minimum bias data sample with an integrated luminos- 
ity of 7.1 fib . Cross sections are measured differen- 
tially in terms of Ar) F , the larger of the pseudorapidity 
regions extending to the limits of the ATLAS sensi- 
tivity, at r\ = ±4.9, in which no final state particles 
are produced above a transverse momentum threshold 
p T ut . The measurements span the region < Arj F < 8 
for 200 < p T ut < 800 MeV. At small A-q F , the data 
test the reliability of hadronisation models in describ- 
ing rapidity and transverse momentum fluctuations in 
final state particle production. The measurements at 
larger gap sizes are dominated by contributions from 
the single diffractive dissociation process (pp — > Xp), 
enhanced by double dissociation (pp — > XY) where 
the invariant mass of the lighter of the two dissocia- 
tion systems satisfies My ~ 7 GeV. The resulting cross 
section is da/dAr/ F w 1 mb for Arf > 3. The large 
rapidity gap data are used to constrain the value of the 
Pomeron intercept appropriate to triple Regge models 
of soft diffraction. The cross section integrated over all 
gap sizes is compared with other LHC inelastic cross 
section measurements. 

PACS 12.40.Nn ■ 12.38.Lg 



1 Introduction 

When two protons collide inelastically at a centre-of- 
mass energy of 7 TeV in the Large Hadron Collider 
(LHC), typically around six charged particles are pro- 
duced with transverse momentunp > 100 MeV per 



1 In the ATLAS coordinate system, the 2-axis points in 
the direction of the anti-clockwise beam viewed from above. 



unit of pseudorapidity in the central region [TH3]. On 
average, the rapidity difference between neighbouring 
particles is therefore around 0.15 units of rapidity, with 
larger gaps occurring due to statistical fluctuations in 
the hadronisation process. Such random processes lead 
to an exponential suppression with gap size [4|, but 
very large gaps are produced where a i-channel colour 
singlet exchange takes place. This may be due to an 
electroweak exchange, but occurs much more frequently 
through the exchange of strongly interacting states. At 
high energies such processes are termed 'diffractive' and 
are associated with 'Pomeron' exchange [5JE] . 

The total cross section in hadronic scattering exper- 
iments is commonly decomposed into four main com- 
ponents: elastic (pp — > pp in the LHC context), single- 
diffractive dissociation (SD, pp — > Xp or pp — > pX, Fig- 
ure[Th.), double-diffractive dissociation (DD, pp — > XY, 
Figure[TjD) and non-diffractive (ND) contributions. The 
more complex central diffractive configuration (CD, 
pp — > pXp, Figure [lfc), in which final state particles are 
produced in the central region with intact protons on 
both sides, is suppressed relative to the SD process by 
a factor of around 10 at high energies [7]. Together, the 
diffractive channels contribute approximately 25-30% 
of the total inelastic cross section at LHC energies [5]. 
Following measurements at the LHC of the elastic [5] , 
total [TU] and total inelastic J5J[TD] cross sections, this 
article contains the first detailed exploration of diffrac- 
tive dissociation processes. 



Polar angles 6 and transverse momenta pt are measured with 
respect to this axis. The pseudorapidity r) = — In tan(0/2) is a 
good approximation to the rapidity of a particle whose mass is 
negligible compared with its energy and is used here, relative 
to the nominal z = point at the centre of the apparatus, to 
describe regions of the detector. 
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(a) (b) (c) 

Fig. 1: Schematic illustrations of the single-diffractive dissociation (a), double-diffractive dissociation (b) and 
central diffractive (c) processes and the kinematic variables used to describe them. By convention, the mass My 
is always smaller than Mx in the double dissociation case and My = M p in the single dissociation case, M p being 
the proton mass. 



Understanding diffractive processes is important in 
its own right, as they are the dominant contribution 
to high energy quasi-elastic scattering between hadrons 
and, via ideas derived from the optical theorem [TTj . 
are also related to the total cross section. They are of- 
ten interpreted at the parton level in terms of the ex- 
change of pairs of gluons [12,13] and are thus sensitive 
to possible parton saturation effects in the low Bjorken- 
x regime of proton structure [T4HI6] . Diffractive cross 
sections also have relevance to cosmic ray physics [17] 
and may be related to the string theory of gravity [18 . 
At the LHC, diffractive dissociation must be well under- 
stood for a good description of the additional inelastic 
proton-proton interactions (pile-up) which accompany 
most events. It also produces a significant uncertainty 
in approaches to luminosity monitoring which rely on 
measurements of the total, or total inelastic, cross sec- 
tion [19] . 

Diffractive dissociation cross sections have been mea- 
sured previously over a wide range of centre-of-mass 
energies. Early measurements are reviewed in [20H24] . 
SD measurements have been made in pp scattering at 
the SPS [25][26] and the Tevatron [27][l8], and also in 
photoproduction [29, 30] and deep inelastic scattering 
[3IHS3] at HERA. Limited high energy DD EgHHJEl] 
and CD [7,35,36 data are also available. In most cases, 
the momentum transfer is too small to permit an inter- 
pretation in terms of partonic degrees of freedom |37] . 
Instead, phenomenological models such as those based 
on Regge theory have been developed [22[l381[3"9"j , which 
underlie the Monte Carlo generators typically used to 
predict the properties of soft inelastic collisions [40H42] . 
Mixed approaches have also been developed which em- 
ploy perturbative QCD where possible [131131]. Large 
theoretical uncertainties remain in the detailed dynam- 
ics expected at the LHC. 



Direct measurements of the masses Mx and My of 
the dissociated systems are difficult at ATLAS, since 
many of the final state particles are produced beyond 
the acceptance of the detector. However, the dissocia- 
tion masses are closely correlated with the size of the ra- 
pidity region in which particle production is completely 
suppressed due to the net colour-singlet Pomeron ex- 
change. This correlation is exploited in this paper, with 
cross sections reported as a function of the size of a 
pseudorapidity region which is devoid of final state par- 
ticle production. These unpopulated pseudorapidity re- 
gions are referred to in the following as 'rapidity gaps', 
or simply 'gaps'. 

To maximise the pseudorapidity coverage and sensi- 
tivity to charged and neutral particle production, rapid- 
ity gaps are identified using both the ATLAS calorime- 
ters and tracking detectors. The specific observable stud- 
ied is Arj F , the larger of the two 'forward' pseudora- 
pidity regions extending to at least 77 = ±4.9 in which 
no particles are produced with pt > jfijP, where p™* 
is varied between 200 MeV and 800 MeV. ND contri- 
butions appear at small gap sizes, with pS 1 * and Arj F 
dependences which are sensitive to fluctuations in the 
hadronisation process. For small p™* choices, the large 
gap size region is dominated by SD events and DD 
events in which one of the dissociation masses is small. 



2 Experimental Method 

2.1 The ATLAS detector 

The ATLAS detector is described in detail elsewhere 
|45j . The beam-line is surrounded by the 'inner detec- 
tor' tracking system, which covers the pseudorapidity 
range |ry| < 2.5. This detector consists of silicon pixel, 
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silicon strip and straw tube detectors and is enclosed 
within a uniform 2 T solenoidal magnetic field. 

The calorimeters lie outside the tracking system. A 
highly segmented electromagnetic (EM) liquid argon 
sampling calorimeter covers the range \rj\ < 3.2. The 
EM calorimeter also includes a pre-sampler covering 
\r)\ < 1.8. The hadronic end-cap (HEC, f .5 < \rj\ < 3.2) 
and forward (FCal, 3.1 < |?7| < 4.9) calorimeters also 
use liquid argon technology, with granularity decreas- 
ing with increasing \rj\. Hadronic energy in the central 
region (|?7| < 1.7) is reconstructed in a steel/scintillator- 
tile calorimeter. 

Minimum bias trigger scintillator (MBTS) detectors 
are mounted in front of the end-cap calorimeters on 
both sides of the interaction point and cover the pseu- 
dorapidity range 2.1 < \r/\ < 3.8. The MBTS is divided 
into inner and outer rings, both of which have eight-fold 
segmentation. The MBTS is used to trigger the events 
analysed here. 

In 2010, the luminosity was measured using a Ceren- 
kov light detector which is located 17 m from the inter- 
action point. The luminosity calibration is determined 
through van der Meer beam scans [I91I46] . 

2.2 Event selection and backgrounds 

The data used in this analysis were collected during the 
first LHC run at y/s = 7 TcV in March 2010, when the 
LHC was filled with two bunches per beam, one pair 
colliding at the ATLAS interaction point. The peak in- 
stantaneous luminosity was 1.1 x 10 27 cm -2 s _1 . Events 
were collected from colliding proton bunch crossings in 
which the MBTS trigger recorded one or more inner 
or outer segments above threshold on at least one side 
of ATLAS. After reconstruction, events are required to 
have hits in at least two of the MBTS segments above 
a threshold of 0.15pC. This threshold cut suppresses 
contributions from noise, which are well modelled by a 
Gaussian with 0.02 pC width. No further event selection 
requirements are applied. 

The data sample corresponds to an integrated lu- 
minosity of 7.1 ± 0.2 /ib" 1 and the number of recorded 
events is 422776. The mean number of interactions per 
bunch crossing is below 0.005, which is consistent with 
the approximately 400 events which have multiple re- 
constructed vertices. Pile-up contamination is thus neg- 
ligible. 

The data sample contains a contribution from beam- 
induced background, mainly due to scattering of beam 
protons from residual gas particles inside the detector 
region. This contamination is estimated using events 
collected in unpaired bunches and is subtracted statis- 
tically in each measurement interval. Averaged over the 



full measurement region, it amounts to 0.2% of the sam- 
ple. More complex backgrounds in which beam-induced 
background is overlaid on a physics event are negligi- 
ble. 



2.3 Reconstruction of rapidity gaps 

The analysis of final state activity in the central region 
(\r]\ < 2.5) is based on combined information from inner 
detector tracks and calorimeter modules. In the region 
2.5 < |?/| < 4.9, beyond the acceptance of the inner de- 
tector, calorimeter information alone is used. The track 
selection is as detailed in [I] . Energy deposits from final 
state particles in the calorimeters are identified using a 
topological clustering algorithm [471148] , with a further 
requirement to improve the control over noise contribu- 
tions, as described below. 

The identification of rapidity gap signatures relies 
crucially on the suppression of calorimeter noise con- 
tributions. The root-mean-squared cell energies due to 
noise vary from around 20 MeV in the most central 
region to around 200 MeV for the most forward re- 
gion [l£r. The shapes of the cell noise distributions in 
each calorimeter are well described by Gaussian distri- 
butions of standard deviation cx no i so , with the excep- 
tion of the tile calorimeter, which has extended tails. 
The default clustering algorithm [48] is seeded by cells 
for which the significance of the measured energy, E, is 
5 = E/cr no \sc > 4. However, with this threshold there 
are on average six clusters reconstructed per empty 
event due to fluctuations in the noise distributions. To 
suppress noise contributions to acceptable levels for gap 
finding, clusters of calorimeter energy deposits are thus 
considered only if they contain at least one cell outside 
the tile calorimeter with an energy significance above 
an ^-dependent threshold, 5th- This threshold is de- 
termined separately in pseudorapidity slices of size 0.1 
such that the probability of finding at least one noisy 
cell in each 77-slice has a common value, 1.4 x 10~ 4 . This 
choice optimises the resolution of the reconstructed gap 
sizes with respect to the gaps in the generated final state 
particle distributions according to MC studies. Since 
the number of cells in an 77-slice varies from about 4000 
in the central region to 10 in the outer part of the FCal, 
the cell thresholds vary between 5th = 5.8 in the cen- 
tral region and 5th = 4.8 at the highest |?7| values in 
the FCal. 

The level of understanding of the calorimeter noise 
is illustrated in Figure [2j which shows the distributions 
of the cell significance 5 for each of the liquid argon 
modules. MBTS-triggered data from colliding bunch 
crossings are compared with a Monte Carlo simulation 
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Fig. 2: Cell energy significance, S = E/a no i ae , distribu- 
tions for the EM (a), HEC (b) and FCal (c) calorime- 
ters. Each cell used in the analysis is included for ev- 
ery event, with the normalisation set to a single event. 
MBTS-triggered minimum bias data (points) are com- 
pared with events randomly triggered on empty bunch 
crossings (histograms) and with a Monte Carlo simula- 
tion (shaded areas). 



and with events which are required to exhibit no ac- 
tivity in the non-calorimeter components of the detec- 
tor, triggered randomly on empty bunch crossings. The 
signal from pp collisions is clearly visible in the long 
positive tails, which are well described by the simula- 
tion. The data from the empty bunch crossings show the 
shape of the noise distribution with no influence from 
physics signals. The empty bunch crossing noise distri- 
butions are symmetric around zero and their negative 
sides closely match the negative parts of the MBTS- 
triggered data distributions. The noise distribution is 
well described over seven orders of magnitude by the 
MC simulation, the small residual differences at pos- 
itive significances being attributable to deficiencies in 
the modelling of pp collision processes. 

The measured energies of calorimeter clusters which 
pass the noise requirements are discriminated using a 
given value of p™S neglecting particle masses. 
The calorimeter energy scale for electromagnetic show- 
ers is determined from electron test-beam studies and 
Z —> e + e~ data [50], confirmed at the relatively small 
energies relevant to the gap finding algorithm through 
a dedicated study of ir° — > 77 decays. The calorimeter 
response to hadronic showers is substantially lower than 
that to electromagnetic showers. In the central region, 
the scale of the hadronic energy measurements is de- 
termined relative to the electromagnetic scale through 
comparisons between the calorimeter and inner detec- 
tor measurements of single isolated hadrons [5TM53) . 
Beyond the acceptance region of the tracking detec- 
tors, the difference between the electromagnetic and 
the hadronic response is determined from test-beam re- 
sults |54H56j . For the purposes of discriminating against 
thresholds in the gap finding algorithm, all cluster en- 
ergy measurements are taken at this hadronic scale. An 
interval in 77 is deemed to contain final state particles 
if at least one cluster in that interval passes the noise 
suppression requirements and has a transverse momen- 
tum above p™*, or if there is at least one good inner 
detector track with transverse momentum above pj? . 

2.4 Definition of forward rapidity gap observable 

The reconstructed forward gap size Ar/ F is defined by 
the larger of the two empty pseudorapidity regions ex- 
tending between the edges of the detector acceptance at 
77 = 4.9 or -q = —4.9 and the nearest track or calorime- 
ter cluster passing the selection requirements at smaller 
|?7|. No requirements are placed on particle production 
at \r)\ > 4.9 and no attempt is made to identify gaps 
in the central region of the detector. The rapidity gap 
size relative to 77 = ±4.9 lies in the range < Ajf < 8, 
such that for example Ar/ F = 8 implies that there is 
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no reconstructed particle with > P™* m one °f the 
regions —4.9 < i] < 3.1 or —3.1 < rj < 4.9. The upper 
limit on the gap size is constrained via the requirement 
of a high trigger efficiency by the acceptance of the 
MBTS detector. 

The measurement is performed in Arj F intervals of 
0.2, except at the smallest values Arj F < 2.0, where 
the differential cross section varies fastest with Ar/ F 
and the gap end-point determination is most strongly 
dependent on the relatively coarse cell granularity of 
the FCal. The bin sizes in this region are increased to 
0.4 pseudorapidity units, commensurate with the reso- 
lution. 

The default value of the transverse momentum thre- 
shold is chosen to be p T ut — 200 MeV. This value lies 
within the acceptance of the track reconstruction for 
the inner detector and ensures that the efficiency of 
the calorimeter cluster selection is greater than 50% 
throughout the r\ region which lies beyond the tracking 
acceptance. 

As described in Section I3T41 the data are fully cor- 
rected for experimental effects using the Monte Carlo 
simulations introduced in Section l3~2l The rapidity gap 
observable defining the measured differential cross sec- 
tions are thus specified in terms of stable (proper life- 
time > 10 ps) final state particles (hereafter referred to 
as the 'hadron level'), with transverse momentum larger 
than the threshold, p T ut , used in the gap reconstruction 
algorithm. 



3 Theoretical Models and Simulations 

3.1 Kinematic Variables and Theory 

As illustrated in Figure [T^, and b, diffractive dissoci- 
ation kinematics can be described in terms of the in- 
variant masses Mx and My of the dissociation sys- 
tems X and y, respectively (with My = M p in the 
SD case), and the squared four-momentum transfer t. 
In the following, the convention My < Mx is adopted. 
The cross section is vastly dominated by small values of 
|i| S; 1 GeV 2 , such that the intact proton in SD events 
is scattered through only a small angle, gaining trans- 
verse momentum px — y/\t\- Further commonly used 
kinematic variables arc defined as 



e 



X 



6 



My 



(1) 



where s is the square of the centre-of-mass energy. 

Diffractive dissociation cross sections can be mod- 
elled using Regge phenomenology pt8l39l!57) , with Pom- 
eron exchange being the dominant process at small £x 
values. For the SD case, the amplitude is factorised into 



a Pomeron flux associated with the proton which re- 
mains intact, and a total probability for the interac- 
tion of the Pomeron with the dissociating proton. The 
latter can be described in terms of a further Pomeron 
exchange using Muller's generalisation of the optical 
theorem [TT], which is applicable for s M x 3> m 2 . 
The SD cross section can then be expressed as a triple 
Pomeron (PPP) amplitude, 

da 



G 31P (0)s 2a ^- 2 (Ml) 



2 \ aip(0)-2a n >(t) 



/(*) (2) 



dtdM x 

where G3ip(0) is a product of couplings and ajp(t) = 
OiTp(Q) + a'jpt is the Pomeron trajectory. The term f(t) 
is usually taken to be exponential such that da/dt oc 
e B(s,M x )t at g xe( j s anc [ Mx, B being the slope pa- 
rameter. With aip(0) close to unity and |t| small, equa- 
tion @ leads to an approximately constant dcr/dln£x 
at fixed s. The DD cross section follows a similar depen- 
dence at fixed £y. The deviations from this behaviour 
are sensitive to the intercept aip(0) of the Pomeron tra- 
jectory |58U59j and to absorptive corrections associated 
with unitarity constraints [431144] . 

The rapidity gap size and its location are closely 
correlated with the variables £x and £y. For the SD 
process, the size Ar\ of the rapidity gap between the 
final state proton and the X system satisfies 

A V ~-]n£ x . (3) 

The Arj F observable studied here differs from Arj in 
that Arj F takes no account of particle production at 
\r\\ > 4.9. For the SD process, where the intact proton 
has rj ~ ±i ln(s/rnp) ~ ±8.9, the gap variables are re- 
lated by Arf ~ Ar\ — 4. Equations © and © thus 
lead to approximately constant predicted cross sections 
dajdArf for SD and low My DD events. With the 
high centre-of-mass energy of the LHC and the exten- 
sive acceptance of the ATLAS detector, events with £x 
between around 10 -6 and 10 -2 can be selected on the 
basis of their rapidity gap signatures, corresponding ap- 
proximately to 7 < M x < 700 GeV. 

Previous proton-proton scattering [25] and photo- 
production [291130] experiments have observed enhance- 
ments relative to triple- Pomeron behaviour at the small- 
est Mx values in the triple Regge region. This effect 
has been interpreted in terms of a further triple Regge 
term (PPR) in which the reaction still proceeds via 
Pomeron exchange, but where the total Pomeron-proton 
cross section is described by a sub-leading Reggeon (JR) 
with intercept «]r(0) — 0.5 [SS]- This leads by analogy 
with equation ([2]) to a contribution to the cross section 
which falls as da/dM x oc 1/M X . In the recent model 
of Ryskin, Martin and Khoze (RMK) g3], a modified 
triple-Pomeron approach to the large £,x region is com- 
bined with a dedicated treatment of low mass diffrac- 
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tive dissociation, motivated by the original s-channel 
picture of Good and Walker [5U], in which proton and 
excited proton eigenstates scatter elastically from the 
target with different absorption coefficients. This leads 
to a considerable enhancement in the low £x cross sec- 
tion which is compatible with that observed in the pre- 
LHC data J25]. 

3.2 Monte Carlo Simulations 

Triple Pomeron-based parameterisations are implemen- 
ted in the commonly used Monte Carlo (MC) event 
generators, PYTHIA [401141] and PHOJET [42]|6l]. These 
generators are used to correct the data for experimental 
effects and as a means of comparing the corrected data 
with theoretical models. 

By default, the PYTHIA model of diffractive dissoci- 
ation processes uses the Schuler and Sjostrand param- 
eterisation [6 2) of the Pomeron flux, which assumes a 
Pomeron intercept of unity and an exponential t depen- 
dence e B ^ x '^ t . Three alternative flux models are also 
implemented. The Bruni and Ingelman version |63) is 
similar to Schuler and Sjostrand, except that its t de- 
pendence is given by the sum of two exponentials. In the 
Berger and Streng [64.65 and Donnachie and Landshoff 
|66] models, the Pomeron trajectory is linear, with vari- 
able parameters, the default being ajp(t) = 1. 085+0. 2bt 
[57], consistent with results from fits to total [551155] 
and elastic [55] hadronic cross section data. Whilst the 
model attributed to Berger and Streng has an exponen- 
tial t dependence, the Donnachie and Landshoff version 
is based on a dipole model of the proton elastic form 
factor. For all flux parameterisations in pythia, addi- 
tional factors are applied to modify the distributions in 
kinematic regions in which a triple-Pomeron approach 
is known to be inappropriate. Their main effects are to 
enhance the low mass components of the dissociation 
spectra, to suppress the production of very large masses 
and, in the DD case, to reduce the probability of the 
systems X and Y overlapping in rapidity space [4T1I62] . 

Above the very low mass resonance region, dissoci- 
ation systems are treated in the PYTHIA6 generator us- 
ing the Lund string model [5J5] , with final state hadrons 
distributed in a longitudinal phase space with limited 
transverse momentum. In PYTHIA8, diffractive parton 
distribution functions from HERA [31] are used to in- 
clude diffractive final states which are characteristic of 
hard partonic collisions, whilst preserving the £xi £yj 
s and t dependences of the diffractive cross sections 
from the PYTHIA6 model [70]. This approach yields a 
significantly harder final state particle transverse mo- 
mentum spectrum in SD and DD processes in PYTHIA8 
compared with PYTHIA6, in better agreement with the 



present data. The default pythia multiple parton in- 
teraction model is applied to ND events and, in the case 
of PYTHIA8, also within the dissociated systems in SD 
and DD events. 

The specific versions used to correct the data are 
PYTHIA6.4.21 (with the AMBT1 tune performed by 
ATLAS [71]) and PYTHIA8.145 (with the 4C tune [72]). 
Updated versions, PYTHIA8.150 and PYTHIA6.4.25 (using 
the 4C and AMBT2B tunes, respectively), are used for 
comparisons with the corrected data (see Table [1]). The 
4C tune of PYTHIA8 takes account of the measurement 
of the diffractive fraction fo of the inelastic cross sec- 
tion in [8], whilst keeping the total cross section fixed, 
resulting in a somewhat smaller diffractive cross section 
than in PYTHIA6. 

The PHOJET model uses the two component dual 
parton model T3J to combine features of Regge phe- 
nomenology with AGK cutting rules [7?] and leading 
order QCD. Diffractive dissociation is described in a 
two-channel eikonal model, combining a triple Regge 
approach to soft processes with lowest order QCD for 
processes with parton scattering transverse momenta 
above 3 GeV. The Pomeron intercept is taken to be 
otjp(Q) — 1.08 and for hard diffraction, the diffractive 
parton densities are taken from [751176] . Hadronisation 
follows the Lund string model, as for pythia. The CD 
process is included at the level of 1.7% of the total in- 
elastic cross section. The specific version used is pho- 
jeti.12.1.35, with fragmentation and hadronisation as 
in PYTHIA6.1.15. 

After integration over t, £x and £y, the cross sec- 
tions for the diffractive processes vary considerably be- 
tween the default MC models, as shown in Table[T] The 
DD variation is particularly large, due to the lack of ex- 
perimental constraints. For use in the data correction 
procedure, the overall fractional non-diffractive (/jvd) 
and diffractive (f D = f SD + f DD + f C D = 1 - Ind) con- 
tributions to the total inelastic cross section are mod- 
ified to match the results obtained in the context of 
each model in a previous ATLAS analysis [8j. Despite 
the close agreement between the diffractive fractions 
fo ~ 30% determined for the three default models (see 
the 'Tuned' fractions in Table [1]), the fu parameter is 
rather sensitive to the choice of Pomeron flux model and 
to the value of aip(0), for example reaching ~ 25% 
for the Bruni and Ingelman flux in PYTHIA8 [5J. 

The default PHOJET and pythia models do not take 
into account Tevatron data which are relevant to the 
decomposition of the diffractive cross section into SD, 
DD and CD components, so these fractions are also ad- 
justed for the present analysis. Based on CDF SD [55] 
and DD [33] cross section data, extrapolated to the 
full diffractive kinematic ranges in each of the mod- 
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Table 1: Predicted ND, SD, DD and CD cross sections, 
together with the fractions of the total inelastic cross 
section Jnd, fsD, Jdd and fcD attributed to each pro- 
cess according to the default versions of the MC models 
(PYTHIA8.150, PYTHIA6.4.25 and PHOJET1.12.1.35), used 
for comparisons with the measured cross sections. The 
modified fractions used in the trigger efficiency and mi- 
gration unfolding procedure, tuned as explained in the 
text, are also given. 



Cross section at y's 


= 7 TeV 




Process 


PYTHIA6 


PYTHIAS 


PHOJET 


ctjvd (mb) 


48.5 


50.9 


61.6 


a SD (mb) 


13.7 


12.4 


10.7 


ctdd (mb) 


9.2 


8.1 


3.9 


<t C d (mb) 


0.0 


0.0 


1.3 


Default f ND (%) 


67.9 


71.3 


79.4 


Default fso (%) 


19.2 


17.3 


13.8 


Default f DD (%) 


12.9 


11.4 


5.1 


Default fen (%) 


0.0 


0.0 


1.7 


Tuned f ND (%) 


70.0 


70.2 


70.2 


Tuned fso (%) 


20.7 


20.6 


16.1 


Tuned fnn (%) 


9.3 


9.2 


11.2 


Tuned fen (%) 


0.0 


0.0 
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els, constraints of 0.29 < ctdd/o'sd < 0.68 and 0.44 < 
&dd/<?sd < 0.94 are derived for the pythia and pho- 
JET models of diffraction, respectively. The tuned ratios 
used in the correction procedure are taken at the cen- 
tres of these bounds. The CD contribution in phojet is 
compatible with the measured Tevatron value of 9.3% 
of the SD cross Section [7] and (Jcd/<?sd is therefore 
kept fixed, with fcD increasing in proportion to Jsd- 
Table [T] summarises the tuned decomposition of the in- 
elastic cross section for each MC model. 

Despite the substantial differences between the ap- 
proaches to diffraction taken in phojet and pythia, 
the two models both employ the Lund string model [BTi] 
of hadronisation. In order to investigate the sensitiv- 
ity of the data at small gap sizes to the hadronisa- 
tion model for ND processes, comparisons of the mea- 
sured cross sections are also made with the HERWIG++ 
generator [77] (version 2.5.1 with the UE7-2 tune [751 
[75]). which uses an alternative cluster-based model. The 
HERWIG++ minimum bias generator takes the total in- 
elastic cross section to be 81 mb, based on a Donnachie- 
Landshoff model [50] • Perturbatively treated semi-hard 
processes are distinguished from soft processes accord- 
ing to whether they produce objects with transverse 
momentum above a fixed threshold which is taken to be 
3.36 GeV. Partons produced from the parton shower are 
combined into colour singlet pairs called clusters, which 
can be interpreted as excited hadronic resonances. The 



clusters are then successively split into new clusters 
until they reach the required mass to form hadrons. 
The most recent HERWIG++ versions contain a mech- 
anism to reconnect partons between cluster pairs via 
a colour reconnection (CR) algorithm, which improves 
the modelling of charged particle multiplicities in pp 
collisions 81 . Similarly to pythia, herwig++ contains 
an eikonalised underlying event model, which assumes 
that separate scatterings in the same event are inde- 
pendent. At fixed impact parameter, this leads to Pois- 
son distributions for both the number of soft scatters 
and the number of semi-hard processes per event. There 
is thus a small probability for 'empty' events to occur 
with no scatterings of either type. Under these circum- 
stances, particle production occurs only in association 
with the dissociation of the beam protons, in a man- 
ner which is reminiscent of diffractive dissociation pro- 
cesses. 

3.3 Comparisons between Monte Carlo simulations 
and uncorrected data 

For use in the correction procedure, MC events are 
processed through the ATLAS detector simulation pro- 
gram [52] . which is based on GEANT4 [53]. They are 
then subjected to the same reconstruction and analysis 
chain as is used for the data. 

The quality of the MC description of the most im- 
portant distributions for the correction procedure is 
tested through a set of control plots which compare the 
uncorrected data and MC distributions. These include 
energy flows, track and calorimeter cluster multiplicities 
and transverse momentum distributions, as well as lead- 
ing cell energy significances in different pseudorapidity 
regions. All such distributions are reasonably well de- 
scribed. Examples are shown in Figures[3^,-d, where the 
total multiplicities of calorimeter clusters which pass 
the selection described in Section 12.31 are shown for 
events in four different regions of reconstructed forward 
rapidity gap size. Whilst none of the MC models gives 
a perfect description, particularly at small multiplici- 
ties, the three models tend to bracket the data, with 
PYTHIA6 showing an excess at low multiplicities and 
PYTHIA8 and phojet showing a deficiency in the same 
region. 

A further example control distribution is shown in 
Figure The probability of detecting at least one 
calorimeter cluster passing the noise requirements with 
Pt > Pt"' = 200 MeV in the most central region < 
0.1) is shown as a function of the px of the leading 
track reconstructed in the same r\ region. In cases where 
this track has pt below around 400 MeV, it spirals in 
the solenoidal field outside the acceptance of the EM 
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Fig. 3: Comparisons of uncorrected distributions between data and MC models, (a)-(d) Total calorimeter cluster 
multiplicities Nc for events reconstructed with (a) < Arf < 2, (b) 2 < Arf < 4, (c) 4 < Arf < 6 and (d) 
6 < Arf < 8. (e) Probability of detecting significant calorimeter energy in the most central region \r)\ < 0.1 as a 
function of the highest transverse momentum max(p^ rack ) of the tracks reconstructed in the inner detector in the 
same \rj\ range. The bin at zero corresponds to events where no charged track with > 160 MeV is reconstructed, 
(f) Forward rapidity gap distribution for p™* = 200 MeV. The final bin at Arf = 10 corresponds to cases where 
no reconstructed particles have pr > p™*- 
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calorimeter. The plotted quantity then corresponds to 
the detection probability for neutral particles in the 
vicinity of a track. Good agreement is observed between 
MC and data. 

The shape of the uncorrected Ar\ F distribution for 
p^ ut = 200 MeV is compared between the data and the 
MC models in Figure |3f. The binning reflects that used 
in the final result (Section 12 .4[) except that contribu- 
tions with Arj F > 8, where the trigger efficiency be- 
comes small, are also shown. None of the models con- 
sidered are able to describe the data over the full Arj F 
range, with the largest deviations observed for small 
non-zero gaps in phojet. All of the models give an 
acceptable description of the shape of the distribution 
for large gaps up to the limit of the measurement at 
Arf = 8 and beyond. 

Considering all control plots together, PYTHIA8 pro- 
vides the best description of the shapes of the distribu- 
tions. Hence this generator is chosen to correct the data. 
The deviations from PYTHIA8 of PYTHIA6 and phojet, 
which often lie in opposite directions and tend to en- 
close the data, are used to evaluate the systematic un- 
certainties on the unfolding procedure. 



3.4 Corrections for Experimental Effects 

After the statistical subtraction of the beam-induced 
background in each interval of Arf (Section [272]) , the 
data are corrected for the influence of the limited accep- 
tance and small particle detection inefficiencies of the 
MBTS using the MC simulation. For the ND, SD and 
DD processes, the trigger efficiency is close to 100% for 
Ar] F < 7, dropping to around 80% at Arf = 8. Since 
the topology of CD events sometimes involves hadronic 
activity in the central region of the detector, with gaps 
on either side, a larger fraction fail the trigger require- 
ment, with efficiencies of close to 100% for Arf < 3 
and between 85% and 95% for 3 < Arf < 8. 

The data are corrected for migrations between the 
reconstructed and hadron level Arf values, due to 
missed or spurious activity and cases where a final state 
particle is observed in a different r\ interval from that 
in which it is produced. The migration corrections are 
obtained using a Bayesian unfolding method [53] with 
a single iteration. The priors for the unfolding proce- 
dure with each MC model are taken after tuning the 
diffractive cross sections as described in Section 13.21 
The migration matrix between the reconstructed and 
hadron level forward gap distributions according to the 
PYTHIA8 MC is shown for p^ ut = 200 MeV in Figure [U 
An approximately diagonal matrix is obtained. 




Fig. 4: Migration matrix between the reconstructed and 
hadron level values of Arf for p™* = 200 MeV, accord- 
ing to PYTHIA8. The distribution is normalised to unity 
in columns and is shown to beyond the limit of the 
measurement at Arj F = 8. 

4 Systematic Uncertainties 

The sources of systematic uncertainty on the measure- 
ment are outlined below. 

MC Model and Unfolding Method Dependence: The trig- 
ger efficiency and migration correction procedure is car- 
ried out using each of the PYTHIA6, PYTHIA8 and pho- 
jet models. The deviation of the data unfolded with 
PHOJET from those obtained with PYTHIA8 is used to 
obtain a systematic uncertainty due to the assumed 
£,x, £,y and t dependences in the unfolding procedure. 
The model dependence due to the details of the final 
state particle production is obtained from the differ- 
ence between the results obtained with PYTHIA6 and 
PYTHIA8. Both of these model dependences are evalu- 
ated separately in each measurement interval and are 
applied symmetrically as upward and downward uncer- 
tainties. They produce the largest uncertainty on the 
measurement over most of the measured range. For 
p™* = 200 MeV, the contributions from the PYTHIA6 
and phojet variations are of similar size. Their com- 
bined effect is typically at the 6% level for large Arf , 
growing to 20% for gaps of around 1.5 pseudorapid- 
ity units. At larger p™* values, the PYTHIA6 source 
becomes dominant. The dependence on the unfolding 
technique has also been studied by switching between 
the default Bayesian method [84] , a method using a Sin- 
gular Value Decomposition of the unfolding matrix |85] 
and a simple bin-to-bin method. The resulting varia- 
tions in the measured cross section are always within 
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the systematic uncertainty denned by varying the MC 
model. 

Modelling of Diffractive Contributions: In addition to 
the differences between the Monte Carlo generators, 
additional systematic uncertainties are applied on the 
modelling of the fractional diffractive cross sections. 
The SD and DD cross sections in the PYTHIA8 model are 
each varied to the limits of the constraints from Teva- 
tron data described in Section 13.21 The fraction Jdd is 
enhanced to 11.3% of the total inelastic cross section, 
with fsD reduced to 18.5% to compensate. At the op- 
posite extreme, fsD is enhanced to 23.2% of the cross 
section, with Jdd reduced to 6.6%. These changes result 
in an uncertainty at the 1% level for Arf > 3. A sys- 
tematic uncertainty on the CD cross section is obtained 
by varying the CD and SD cross sections in PHOJET 
between the tuned values and o'cd/o'sd — 0.093, corre- 
sponding to the CDF measurement in [7| . This variation 
also results in a 1% uncertainty in the large gap region. 

Calorimeter Energy Scale: The uncertainty on the calo- 
rimeter energy scale is constrained to be below the 5% 
level down to energies of a few hundred MeV in the 
central region, \n\ < 2.3, through comparisons between 
isolated calorimeter cluster energy measurements and 
momentum determinations of matched tracks in the in- 
ner detector [5TM53] . This method is not available for 
larger \r/\ values beyond the tracking acceptance. How- 
ever, as 1 77 1 grows, the default p™* = 200 MeV thresh- 
old corresponds to increasingly large energies, reaching 
beyond 10 GeV at the outer limits of the FCal. The 
uncertainty on the response to electromagnetic show- 
ers in this energy range is determined as a function of 
1 77 1 from the maximum observed deviations between the 
data and the MC simulation in the peaks of 7r° — > 77 
signals, under a variety of assumptions on background 
shapes and cluster energy resolutions. The relative re- 
sponse to charged pions compared with the electromag- 
netic scale has been studied in the relevant energy range 
for the FCal [55l|86] and HEC [54l[86] test-beam data, 
with systematic uncertainties of 8% and 4%, respec- 
tively, determined from the difference between data and 
MC. Adding the uncertainties in the electromagnetic 
scale and in the relative response to hadrons in quadra- 
ture, energy scale uncertainties of 5% for I77I < 1.37, 
21% for 1.37 < |?7| < 1.52 (transition region between 
barrel and end-cap), 5% for 1.52 < \r)\ < 2.3, 13% for 
2.3 < \r)\ < 3.2 and 12% for 3.2 < \r)\ < 4.9 are ascribed. 
In addition to the absolute calorimeter response, these 
values account for systematic effects arising from dead 
material uncertainties and from the final state decom- 
position into different particle species. In the unfold- 
ing procedure, the corresponding systematic variation 



is applied to energy depositions from simulated final 
state particles in MC, with noise contributions left un- 
changed. The clustering algorithm is then re-run over 
the modified calorimeter cells. The scale uncertainty 
variation is thus considered both in the application of 
the p™* threshold to the clusters and in the discrimina- 
tion of cells within selected clusters against the signif- 
icance cut used to veto noise. The resulting fractional 
uncertainties on the differential cross sections at the de- 
fault p| ut = 200 MeV are largest (reaching ~ 12%) in 
the region Arj F ^ 3, where the gap identification re- 
lies most strongly on the calorimeter information. For 
larger gaps, the well measured tracks play an increas- 
ingly important role in defining the gap size and the 
cross section is dominated by low diffractive events 
for which particle production in the gap region is com- 
pletely suppressed. The sensitivity to the calorimeter 
scale is correspondingly reduced to a few percent. 

MBTS Efficiency: The description of the MBTS effi- 
ciency in the MC models leads to a potential system- 
atic effect on the trigger efficiency and on the off-line 
MBTS requirement. Following [5], the associated uncer- 
tainty is evaluated by increasing the thresholds of all 
MBTS counters in the simulation to match the max- 
imum variation in the measured response in data ac- 
cording to studies with particles extrapolated from the 
tracker or FCal. This systematic error amounts to typ- 
ically 0.5 — 1% for Ar] F > 2 and is negligible at the 
smallest Ar/ F . 

Tracking Efficiency: The dominant uncertainty in the 
charged particle track reconstruction efficiency arises 
due to possible inadequacies in the modelling of the 
material through which the charged particles pass pQ. 
This uncertainty is quantified by studying the influence 
on the data correction procedure of using an MC sample 
produced with a 10% enhancement in the support ma- 
terial in the inner detector. The resulting uncertainty 
is smaller than 3.5% throughout the measured distri- 
bution. 

Luminosity: Following the van der Meer scan results 
in [46j , the normalisation uncertainty due to the uncer- 
tainty on the integrated luminosity is 3.4%. 

Each of the systematic uncertainties is determined 
with correlations between bins taken into account in the 
unfolding by repeating the full analysis using data or 
MC distributions after application of the relevant sys- 
tematic shift. The final systematic error on the differen- 
tial cross section is taken to be the sum in quadrature 
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of all sources. Compared with the systematic uncertain- 
ties, the statistical errors are negligible at the smallest 
gap sizes, where the differential cross section is largest. 
For gap sizes Arj F > 3, the statistical errors are at the 
1% level and are typically smaller than the systematic 
errors by factors between five and ten. 



5 Results 

5.1 Differential cross section for forward rapidity gaps 

In this section, measurements are presented of the in- 
elastic cross section differential in forward rapidity gap 
size, Ar\ F , as defined in Section HOI The data cover the 
range < Arj F < 8. In the large gap region which is 
populated by diffractive processes, the cross section cor- 
responds to a ^-integrated sum of SD events in which ei- 
ther of the colliding protons dissociates and DD events 
with £_ Y Z 1CT 6 {M Y S 7 GeV). The data span the 
range £x ~ 10~ 5 . Diffractive events with smaller £x 
values are subject to large MBTS trigger inefficiencies 
and thus lie beyond the kinematic range of the mea- 
surement. 

As discussed in Section 12.41 the lowest transverse 
momentum requirement for the gap definition which is 
directly accessible experimentally is p™* = 200 MeV. 
Figure [SJi shows the differential gap cross section for 
this choice of p™ , which is also given numerically in 
Table [5J The uncertainty on the measurement is typ- 
ically less than 8% for Arf > 3, growing to around 
20% at Ai] F = 1.5 before improving to around 10% for 
events with little or no forward gap. The data are com- 
pared with the predictions of the default settings of the 
pythia6 (labelled 'pythia6 atlas ambt2b') pythias 
('pythiaS 4C') and phojet models. In Figures [5b-d, 
the results are compared with each of the MC models 
separately, with the default decomposition of the cross 
section into ND, SD, DD and CD contributions accord- 
ing to the models (Table [1} also indicated. 



5.2 Small gap sizes and constraints on hadronisation 
models 

At Arf < 2, all models agree that the ND process is 
dominant and the expected [3] exponential decrease of 
the cross section with increasing gap size, characteristic 
of hadronisation fluctuations, is the dominant feature of 
the data. According to the models, this region also con- 
tains DD events which have £y S 10~ 6 , such that the 
Y system extends into the ATLAS detector acceptance, 
as well as both SD and DD events with very large £x, 



such that no large rapidity gap is present within the re- 
gion \tj\ < 4.9. The default MC models tend to lie above 
the data in this region, a result which is consistent with 
the overestimates of the total inelastic cross section ob- 
served for the same models in [8]. The PYTHIA8 model 
is closest in shape to the data, which is partly due to 
the modification of fr> in the most recent versions made 
in light of the previous ATLAS data [5] • Both PYTHIA 
models are closer to the small Arf data than phojet, 
which exhibits an excess of almost a factor of two for 
Ar) F ~ 1. 

As can be inferred from comparisons between the 
predicted shapes of the ND contributions in the dif- 
ferent MC models (Figures EJ>d), there are consider- 
able uncertainties in the probability of obtaining large 
hadronisation fluctuations among low transverse mo- 
mentum final state particles [8 7) . Studying the depen- 
dence of the measured differential cross section on 
provides a detailed probe of fluctuations in the hadro- 
nisation process in soft scattering and of hadronisation 
models in general. The measurement is thus repeated 
with different choices of pjP , applied both in the rapid- 
ity gap reconstruction and in the definition of the mea- 
sured hadron level cross section. To avoid cases where 
the largest gap switches from one side of the detec- 
tor to the other when low pt particles are excluded by 
the increased p™* choice, the side of the detector on 
which the gap is located is fixed to that determined at 
p™* = 200 MeV for all measured cross sections. 

A comparison between the results with p^S 1 ' — 200, 
400, 600 and 800 MeV is shown in FigureHi. Figures[6p- 
d show the results for p c ^ lt = 400, 600 and 800 MeV, 
respectively, compared with the PYTHIA8, PYTHIA6 and 
phojet MC models. The ND contributions according 
to each of the models are also shown. As p^ ut increases, 
the exponential fall becomes less steep, so larger Arf 
values become more heavily populated and the non- 
diffractive and diffractive contributions in the models 
become similar. Also, the uncertainties due to the MC 
model dependence of the unfolding procedure grow. 

The influence of changing from p^ — 200 MeV 
to pj?* — 400 MeV is small at large Arf , where the 
cross section is dominated by small £x diffractive events 
and particle production is kinematically forbidden over 
a wide range of pseudorapidity. For Ai] F < 4, where 
ND contributions become important, a significant frac- 
tion of events are assessed as having larger gaps for 
p™* = 400 MeV than for p T ut = 200 MeV. As the value 
of p| ut increases to 600 and 800 MeV, soft ND events 
migrate to larger Ai] F values, giving significant contri- 
butions throughout most of the distribution and con- 
firming [T] that the production of final state particles 
with more than a few hundred MeV is rare in minimum 
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Fig. 5: Inelastic cross section differential in forward gap size Arj F for particles with pt > 200 MeV. The shaded 
bands represent the total uncertainties. The full lines show the predictions of PHOJET and the default versions of 
PYTHIA6 and PYTHIA8 The dashed lines in (b-d) represent the contributions of the ND, SD and DD components 
according to the models. The CD contribution according to PHOJET is also shown in (d). 



bias events, even at LHC energies. All MC models are 
able to reproduce the general trends of the data, though 
none provides a full description. 

It is interesting to investigate the extent to which 
the alternative cluster-based approach to hadronisation 
in the non-diffractive HERWIG++ model is able to de- 
scribe the data at small gap sizes, where the contribu- 
tion from ND processes is dominant. A comparison of 
the data at each of the p^ ut values with herwig++ is 
shown in Figure Four versions of the UE7-2 tune are 
shown, with variations in the details of the model which 
are expected to have the largest influence on rapidity 
gap distributions. These are the default version (UE7- 
2), a version in which the colour reconnection model 



is switched off (UE7-2, No CR) and similar versions 
which exclude events with no scatterings of either the 
soft or semi-hard types (UE7-2, No Empty Evts and 
UE7-2, No Empty Evts, No CR). At small gap sizes, all 
versions of the model produce an exponential fall with 
increasing gap size, though the dependence on Arf is 
not steep enough in the default model and is too steep 
when colour recombination effects are switched off. 

Despite not containing an explicit diffractive com- 
ponent, the default HERWIG++ minimum bias model 
produces a sizeable fraction of events with large gaps, 
overshooting the measured cross section by up to factor 
of four in the interval 2 < Arj F < 7 and producing an 
enhancement centred around Arf = 6. When colour 
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Fig. 6: Inelastic cross section differential in forward gap size Ar\ F for different p™' values, (a) Comparison between 
the measured cross sections. The full uncertainties are shown. They are correlated between the different p^ lt choices, 
(b-d) Comparison between the data and the MC models for p™* = 400, 600 and 800 MeV. The non-diffractive 
component in each MC model is also shown. 



reconnection is switched off, this large gap contribution 
is reduced considerably, but remains at a similar level to 
that measured in the range 3 < Arj F < 5. The enhance- 
ment near Ar] F 6 is still present. The events with 
zero scatters in the HERWIG++ underlying event model 
provide a partial explanation for the large gap contribu- 
tion. Removing this contribution reduces the predicted 
large gap cross section, but the non-exponential tail 
and large Arf enhancement persist. For all scenarios 
considered, the alternative cluster based hadronisation 
model in HERWIG++ shows structure which is incom- 
patible with the data. 



5.3 Large gap sizes and sensitivity to diffractive 
dynamics 

At large Arf , the differential cross section exhibits a 
plateau, which is attributed mainly to diffractive pro- 
cesses (SD events, together with DD events at £y ^ 
10~ 6 ) and is shown in detail in Figure [5] According 
to the phojet MC model, the CD contribution is also 
distributed fairly uniformly across this region. Over a 
wide range of gap sizes with Ar/ F > 3, the differential 
cross section is roughly constant at around 1 mb per 
unit of rapidity gap size. Given the close correlation 
between Arf and — In £ (Section 13. ip , this behaviour 
is expected as a consequence of the dominance of soft 
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Fig. 7: Inelastic cross section differential in forward gap size Arj F for pJp* = (a) 200 MeV, (b) 400 MeV, (c) 
600 MeV and (d) 800 MeV. The data are compared with the UE7-2 tune of the HERWIG++ model. In addition to 
the default tune, versions are shown in which the colour reconnection model is switched off and in which events 
with zero scatters are excluded (see text for further details). 



diffractive processes. All MC models roughly reproduce 
the diffractive plateau, though none gives a detailed de- 
scription of the shape as a function of Ai] F . 

When absolutely normalised, the PYTHIA predic- 
tions overshoot the data throughout most of the diffrac- 
tive region, despite the tuning of fjj to previous ATLAS 
data [5] in these models. The excess here is partially a 
reflection of the 10% overestimate of the pythia pre- 
diction in the total inelastic cross section and may also 
be associated with the large DD cross section in the 
measured region, which exceeds that expected based 
on Tevatron data [M] and gives rise to almost equal SD 
and DD contributions at large Ar/ F . For phojet, the 
underestimate of the diffractive fraction fr> is largely 



compensated by the excess in the total inelastic cross 
section, such that the large gap cross section is in fair 
agreement with the measurement up to Arj F w 6. The 
DD contribution to the cross section in phojet is heav- 
ily suppressed compared with that in the pythia mod- 
els. 

Integrated over the diffractive-dominated region 
5 < Ar\ F < 8, corresponding approximately to 
—5.1 ;$ log 10 ;$ —3.8 according to the MC models, 
the measured cross section is 3.05 ± 0.23 mb, approxi- 
mately 4% of the total inelastic cross section. This can 
be compared with 3.58 mb, 3.89 mb and 2.71 mb for 
the default versions of PYTHIA8, PYTHIA6 and phojet, 
respectively. 
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Fig. 8: Inelastic cross section differential in forward gap size Arj F for particles with px > 200 MeV and Ar/ F > 2. 
The error bars indicate the total uncertainties. In (a), the full lines show the predictions of phojet, the default 
versions of PYTHIA6 and PYTHIA8, and PYTHIA8 with the Donnachie-Landshoff Pomcron flux. The remaining plots 
show the contributions of the SD, DD and ND components according to each generator. The CD contribution 
according to phojet is also shown in (d). 



As can be seen in Figure [5J the differential cross 
section rises slowly with increasing Ar/ F for Arj F > 5. 
Non-diffractive contributions in this region are small 
and fall with increasing Ar/ F according to all models, 
so this rise is attributable to the dynamics of the SD 
and DD processes. Specifically the rising cross section is 
as expected from the PPP term in triple Regge mod- 
els with a Pomeron intercept in excess of unity (see 
equation @). In Figure^, a comparison is made with 
the PYTHIA8 model, after replacing the default Schuler 
and Sjostrand Pomeron flux with the Donnachie and 
Landshoff (DL) version using the default Pomeron tra- 
jectory, ajp{t) = 1.085+0.25* ('pythias dl'). It is clear 



that the data at large Arj F are not perfectly described 
with this choice. 

Whilst the data are insensitive to the choice of ap/, 
there is considerable sensitivity to the value of ajp(Q). 
The data in the cleanest diffractive region Ar/ F > 6 
are used to obtain a best estimate of the appropriate 
choice of the Pomeron intercept to describe the data. 
SD and DD PYTHIA8 samples are generated with the 
DL Pomeron flux for a range of ap(0) values. In each 
case, the default a^pi value of 0.25 GeV - is taken and 
the tuned ratios of the SD and DD contributions appro- 
priate to PYTHIA8 from Table Q] are used. The x 2 value 
for the best fit to the data in the region 6 < Ar] F < 8 
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is obtained for each of the samples with different ap(0) 
values, with the cross section integrated over the fitted 
region allowed to float as a free parameter. The opti- 
mum a]p(0) is determined from the minimum of the 
resulting \ 2 parabola. 

The full procedure is repeated for data points shifted 
according to each of the systematic effects described in 
SectionQJ such that correlations between the uncertain- 
ties on the data points are taken into account in eval- 
uating the uncertainties. The systematic uncertainty is 
dominated by the MC model dependence of the data 
correction procedure, in particular the effect of unfold- 
ing using PYTHIA6 in place of PYTHIA8, which leads to 
a significantly flatter dependence of the data on Arf 
at large gap sizes. 

The result obtained in the context of the PYTHIA8 
model with the DL flux parameterisation is 



aip(O) = 1.058 ± 0.003(stat.)tao3g(syst.) 



(4) 



The data are thus compatible with a value of ap(0) 
which matches that appropriate to the description of 
total hadronic cross sections [551155] . When the Berger- 
Streng Pomeron flux, which differs from the DL ver- 
sion in the modelling of the t dependence, is used in 
the fit procedure, the result is modified to ap(0) = 
1.056. The effects of varying aw' between 0.1 GeV -2 
and 0.4 GeV -2 and of varying the /sd and foD frac- 
tions assumed in the fit in the ranges given in Section 2] 
are also smaller than the statistical uncertainty. Com- 
patible results are obtained by fitting the higher p™* 
data. 

A comparison between the data and a modified ver- 
sion of PYTHIA8, with a]p(0) as obtained from the fit, 
is shown in Figure [5] Here, the diffractive contribution 
to the inelastic cross section fn = 25.6% is matched^) 
to the fitted value of ap (0) using the results in [5] . To- 
gether with the cross section integrated over the region 
6 < Aif < 8 as obtained from the fit and the tuned ra- 
tio foD I fsD from Table[U this fixes the normalisation of 
the full distribution. The description of the data at large 
Arf is excellent and the exponential fall at small Ar\ F 
is also adequately described. There is a discrepancy in 
the region 2 < Ar\ F < 4, which may be a consequence 
of the uncertainty in modelling large hadronisation fluc- 
tuations in ND events (compare the ND tails to large 
Arf in Figured, c and d). It may also be attributable 
to sub-leading trajectory exchanges [2HHHD or to the 
lack of a CD component in the pythia model. 

2 Since only data at large Arf are included in the fit, the 
result for ajp (0) is insensitive to systematic variations in fo ■ 
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Fig. 9: Inelastic cross section differential in forward gap 
size Ar] F for particles with pt > 200 MeV. The data 
are compared with a modified version of the PYTHIA8 
model with the DL flux, in which the Pomeron intercept 
aip (0) is determined from fits to the data in the region 
6 < Arf < 8. See text for further details. 

5.4 The integrated inelastic cross section 

By summing over the Arf distribution from zero to 
a maximum gap size Arj F ut , the integrated inelastic 
cross section can be obtained, excluding the contribu- 
tion from events with very large gaps Arf > ArjQ Ut . As 
discussed in Section l3~Tl there is a strong correlation be- 
tween the size of the gap and the kinematics of diffrac- 
tion (see e.g. equation ^ for the SD process). The cross 
section integrated over a given range of gap size can 
thus be converted into an integral over the inelastic pp 
cross section down to some minimum value £cut of £x ■ 
The variation in the integrated inelastic cross section 
with Ar/ F ut can then be used to compare inelastic cross 
section results with different lower limits, £cut- 
The integral of the forward gap cross section 

" AvLt da , A F 
dArj F 



10 W 

is obtained for Ar/ F ut values varying between 3 and 8 
by cumulatively adding the cross section contributions 
from successive bins of the measured gap distribution. 
The correspondence between maximum gap size and 
minimum £x used here is determined from the PYTHIA8 
model to be log 10 £cut — ~0.45Z\^c ut — 1.52. The un- 
certainty on this correlation is small; for example the 
phojet model results in the same slope of —0.45 with 
an intercept of —1.56. This correlation is applied to 
convert to an integral 

/ -JF~ d ^ x ■ 
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A small correction is applied to account for the fact 
that the gap cross section neglects particles witfjfl px < 
p<j, ut = 200 MeV and includes a contribution from ND 
processes. This correction factor is calculated using 
PYTHIA8 with the DL flux, and the optimised ap(0) 
and /d values, as determined in Section 15.31 The inte- 
gration range is chosen such that the correction is al- 
ways smaller than ±1.3%. The systematic uncertainty 
on the correction factor, evaluated by comparison with 
results obtained using phojet or PYTHIA8 with the de- 
fault Schuler and Sjostrand flux, together with the sys- 
tematic variations of the tuned fractions fsD and Jdd 
as in Section |H is also small. 

The integrated inelastic cross section is shown as a 
function of £cut in Figure [TUI where it is also compared 
with a previous ATLAS result [SJ and with the TOTEM 
extraction of the full inelastic cross section [TU], de- 
rived from a measurement of the elastic cross section 
via the optical theorem. The errors on all of the ex- 
perimental data points are dominated by the luminos- 
ity uncertainties. The previous ATLAS result was also 
based on MBTS-triggered data, but is quoted at the 
£x value corresponding to 50% trigger efficiency, which 
is slightly beyond the range accessed here. Extrapolat- 
ing according to the measured dependence on £cut, the 
new data are in good agreement with the previous re- 
sult, the small apparent difference being well within the 
uncertainty due to run-to-run luminosity measurement 
variations. 

It is instructive to compare the TOTEM result with 
the ATLAS measurements, since the latter omit the 
poorly understood lowest £x region. By comparing the 
lowest £cut data point from the present analysis with 
the TOTEM measurement and neglecting any correla- 
tions between the ATLAS and TOTEM uncertainties, 
the inelastic cross section integrated over 
£x < 8 X 10~ 6 is inferred to be 14.5^;° mb. Signif- 
icantly smaller contributions are predicted by the de- 
fault versions of pythia (~ 6 mb) and PHOJET 
(~ 3 mb). FigurcfTUlalso shows two versions of the RMK 
model (see Section 13. ip , corresponding to versions (i) 
(upper curve) and (ii) (lower curve) in |43j . These ver- 
sions differ in the radii attributed to the elastically scat- 
tered eigenstates comprising the low £x contribution 
which is added to the more standard triple Pomeron 

3 The finite pSj, ut value in the measured gap cross sections 
tends to increase gap sizes slightly relative to pij, ut = 0. How- 
ever, MC studies indicate that this effect has the biggest in- 
fluence on the exponentially falling distribution at small gap 
sizes, whereas the difference for the Arj F values which are 
relevant to the integrated cross section are relatively small. 
According to the MC models, the cross section integrated 
over 5 < Ar\ F < 8 decreases by 2% when changing from 
p™* = 200 MeV to p^ ut = 0. 
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Fig. 10: Inelastic cross section excluding diffractive pro- 
cesses with l;x < £cut, obtained by integration of the 
differential cross section from gap sizes of zero to a vari- 
able maximum. The results from the present analysis 
('ATLAS L = 7.1 jitb ) are compared with a previous 
ATLAS result [5J and a TOTEM measurement inte- 
grated over all kinematically accessible £x values [TP] . 
The predictions of the default versions of the PYTHIA6, 
PYTHIA8 and phojet models are also shown, along with 
two versions of the RMK model [33] (see text). The ver- 
tical error bars on the ATLAS measurements denote the 
systematic uncertainties excluding that on the luminos- 
ity measurement, whilst the shaded area represents the 
full systematic uncertainty. For the TOTEM point, the 
error bar represents the statistical uncertainty whilst 
the shaded area represents the full uncertainty. The un- 
certainties on the data points obtained in the present 
analysis are strongly correlated between neighbouring 
points. 



calculation in the model, (ii) being the favoured version 
and (i) being indicative of the flexibility in the model 
whilst preserving an acceptable description of pre-LHC 
data. The additional low t;x processes enhance the in- 
elastic cross section by 5.5 mb and 6.7 mb in versions 
(i) and (ii) , respectively. Although the RMK model lies 
below the data in general, the low £x enhancement is 
compatible with that observed [55]. The shape of the 
distribution at low ^x is not predicted in the model, 
but is compatible with the data if, as shown here |89) . 
it is assumed to have the steep £x dependence associ- 
ated with the PPP, rather than the PPP triple Regge 
term. Similar conclusions have been reached previously 
from proton-proton [35] and photoproduction [2l?Il3"0"] 
data. 
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6 Summary 

A novel algorithm has been devised for identifying ra- 
pidity gaps in the final state of minimum bias ATLAS 
data, leading to measurements in which particle pro- 
duction is considered down to transverse momentum 
thresholds p T ut between 200 MeV and 800 MeV. The 
differential cross section d<r/dAi] F is measured for for- 
ward rapidity gaps of size < Arf < 8, corresponding 
to the larger of the two gaps extending to rj = ±4.9, 
with no requirements on activity at \t}\ > 4.9. An expo- 
nentially falling non-diffractive contribution is observed 
at small gap sizes, which is also a feature of the pythia, 
phojet and HERWIG++ Monte Carlo models. However, 
none of the models describes the Ai] F or p™* depen- 
dence of this region in detail. At large gap sizes, the 
differential cross section exhibits a plateau, which cor- 
responds to a mixture of the single-diffractive dissoci- 
ation process and double dissociation with £y < 10~ 6 . 
This plateau amounts to a cross section close to 1 mb 
per unit of gap size and its magnitude is roughly de- 
scribed by the pythia and phojet Monte Carlo mod- 
els. None of the default models reproduce the rise of the 
differential cross section as a function of gap size at the 
largest Ar\ F values. This rise is interpreted within the 
triple Pomeron-based approach of the PYTHIA8 model 
with a Donnachie-Landshoff Pomeron flux in terms of 
a Pomeron intercept of ap(0) = 1.058 ± 0.003(stat.) 
iao39( s y s t-)- Since the bulk of the inelastic pp cross 
section is contained within the measured range, inte- 
grated cross sections are also obtained and compared 
with previous measurements. The contribution to the 
total inelastic cross section from the region £x < 10~ 5 
is determined to be around 20%, which is considerably 
larger than is predicted by most models. 
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Table 2: The measured differential cross section data points for p™* = 200 MeV, with each value corresponding 
to an average over the given Arf range. Also quoted are the percentage statistical (<$stat) uncertainty and the 
upward (+<$tot) and downward (— <5tot) total uncertainties, obtained by adding all uncertainties in quadrature. 
The remaining columns contain the percentage shifts due to each of the contributing systematic sources, which are 
correlated between data points. Those due to the modelling of final state particle production (^pyg), the modelling 
of the £x, £,y and t dependences (<5p n0 ) and variation of the CD (# C( j) cross section in the unfolding procedure are 
applied symmetrically as upward and downward uncertainties, as are those due to the dead material budget in the 
tracking region (<5 m at) and the MBTS response (<^ m bts) - ^ ne uncertainties due to variations in the relative energy 
scale in data and MC are evaluated separately for upward (<5 e +) and downward (<5 e _) shifts, as are the modelling 
uncertainties due to enhancing (S s a) or reducing ((^h) the <jsd/&dd cross section ratio. Minus signs appear where 
the shift in a variable is anti-correlated rather than correlated with the shift in the differential cross section. The 
3.4% normalisation uncertainty due to the luminosity measurement is also included in the ±<$tot values. These 
data points can be obtained from the HEPDATA database [90], along with their counterparts for p™* = 400, 600 
and 800 MeV. A Rivet 91 routine is also available. 
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